Spectral energy dynamics in magnetohydrodynamic turbulence 
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Spectral direct numerical simulations of incompressible MHD turbulence at a resolution of up to 
1024 3 collocation points are presented for a statistically isotropic system as well as for a setup with 
an imposed strong mean magnetic field. The spectra of residual energy, Ef' = \E% — E^\, and total 
energy, E^ = E^ + E^ 1 , are observed to scale self-similarly in the inertial range as ~ fc _7//3 , 
Ek ~ fc _5//3 (isotropic case) and E^ ~ kj 2 , Ek x ~ fc^ 3//2 (anisotropic case, perpendicular to the 
mean field direction). A model of dynamic equilibrium between kinetic and magnetic energy, based 
on the corresponding evolution equations of the eddy-damped quasi-normal Markovian (EDQNM) 
closure approximation, explains the findings. The assumed interplay of turbulent dynamo and 
Alfven effect yields E^} ~ kE 2 . which is confirmed by the simulations. 



The nonlinear behavior of turbulent plasmas gives rise 
to a variety of dynamical effects such as self-organization 
of magnetic confinement configurations in laboratory ex- 
periments J]|, generation of stellar magnetic fields |2| or 
structure formation in the interstellar medium 0]. The 
understanding of these phenomena is incomplete as the 
same is true for many inherent properties of the under- 
lying turbulence. 

Large-scale low-frequency plasma turbulence is treated 
in the magnetohydrodynamic (MHD) approximation de- 
scribing the medium as a viscous and electrically resis- 
tive magnetofluid neglecting additional kinetic effects. 
Incompressiblity of the flow is assumed for the sake of 
simplicity. In this setting the nature of the turbulent 
energy cascade is a central and still debated issue with 
different phenomenologies being proposed 0, H, IE S 
(cf. for a review). The associated spectral dynamics 
of kinetic and magnetic energy, in spite of its comparable 
importance, has received less attention (as an exception 
see H3). 

This Letter reports a spectral relation between residual 
and total energy, Ef = \Ef - Ef\ and E k = Ef + E^ 
respectively, as well as the influence of an imposed mean 
magnetic field on the spectra. The proposed physical pic- 
ture, which is confirmed by accompanying direct numer- 
ical simulations, embraces two-dimensional MHD tur- 
bulence, globally isotropic three-dimensional systems as 
well as turbulence permeated by a strong mean magnetic 
field. 

In the following reference is made to two high- 
resolution pseudospectral direct numerical simulations 
of incompressible MHD turbulence which we regard as 
paradigms for isotropic (I) and anisotropic (II) MHD tur- 
bulence. The dimensionless MHD equations 
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are solved in a 27r-periodic cube with spherical mode 
truncation to reduce numerical aliasing errors 11]. The 
equations include the flow vorticity, uj = V x v, the mag- 
netic field expressed in Alfven speed units, b, as well as 
dimensionless viscosity, /i, and resistivity, r\. In simula- 
tion II forcing is applied by freezing the largest spatial 
scales of velocity and magnetic field. 

Simulation I evolves globally isotropic freely decaying 
turbulence represented by 1024 3 Fourier modes. The ini- 
tial fields are smooth with random phases and fluctuation 
amplitudes following exp(— k 2 l(2k^f) with fco = 4. To- 
tal kinetic and magnetic energy are initially equal with 
e k _ e m = 5 The ratio e k /E m decreases in time 
taking on values of 0.28 — 0.23 in the period consid- 
ered (cf. [I2]). Tlie- ratio of kinetic and magnetic en- 
ergy dissipation rate, e K je M , with fj, = rj = 1 x 10~ 4 
also decreases during turbulence decay from 0.7 to about 
0.6, the difference in dissipation rates reflecting the im- 
balance of the related energies. The Reynolds number 
Re= (£ K )£y(^£ total ) at t = 6 is about 2700 and slightly 
diminishes during the run. Magnetic, H M = \ J y dVa-b, 
b = V x a, and cross helicity, H c = | J v dVv ■ b , are 
negligible with H c showing a dynamically unimportant 
increase from 0.03 to 0.07 during the simulation. The 
run covers 9 eddy turnover times defined as the time re- 
quired to reach the maximum of dissipation from t = 0. 
The large-scale rms magnetic field decays from initially 
0.7 to 0.3. 

Case II is a 1024 2 x 256 forced turbulence simula- 
tion with an imposed constant mean magnetic field of 
strength bg = 5 in units of the large-scale rms magnetic 
field 6 rms ~ v rms ~ 1. The forcing, which keeps the ratio 
of fluctuations to mean field approximately constant, is 
implemented by freezing modes with k < kf = 2. The 
simulation with fj, = r\ = 9 x 10 -5 has been brought into 
quasi-equilibrium over 20 eddy-turnover times at a res- 
olution of 512 2 x 256 and spans about 5 eddy turnover 
times of quasi-stationary turbulence with 1024 2 x 256 
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FIG. 1: Total (solid), kinetic (dashed), and magnetic (dotted) 
energy spectra in 1024 3 case I simulation (normalized, time- 
averaged and compensated). Dash-dotted line: k~ 3 ^ 2 scaling. 

Fourier modes and Re«2300 (based on field perpendicu- 
lar fluctuations). Kinetic and magnetic energy as well 
as the ratio E K /E M are approximately unity with a 
slight excess of E M . Perpendicular to the imposed field, 
large-scale magnetic fluctuations with & rms ~ 0.4 are ob- 
served. Correspondingly, e K je M ~ 0.95 during the sim- 
ulation. The system has relaxed to H c ~ 0.15 with a 
fluctuation level of about 30% and H M ~ 0.2ff^ ax with 
H Max ~ EM /kf- 

Fourier-space-angle integrated spectra of total, mag- 
netic, and kinetic energy for case I are shown in Fig. ^ 
To neutralize secular changes as a consequence of turbu- 
lence decay, amplitude normalization is used assuming 
a Kolmogorov total energy spectrum, — > Ek/(s(J> 5 ), 
e = —dtE, with wavenumbers given in inverse multi- 
ples of the associated dissipation length, Id ~ (/i 3 /^) 1 / 4 . 
The quasi-stationary normalized spectra are time aver- 
aged over the period of self-similar decay, t = 6 — 8.9. 
As in previous numerical wo rk |l3l|l4ll and also observed 
in solar wind measurements |l5l llfij , Kolmogorov scaling 
applies for the total energy in the well-developed iner- 
tial range, 0.01 < k < 0.1. However, here the remark- 
able growth of excess magnetic energy with decreasing 
wavenumber is of interest. Qualitatively similar behavior 
is observed with large scale forcing exerted on the system. 
We note that no pile-up of energy is seen at the dissipa- 
tive fall-off contrary to other high-resolution simulations 
[Til \v\ . Apart from different numerical techniques and 
physical models this difference might be due to the lim- 
ited simulation period at highest resolution namely 5 [l^j 
and 4.3 large-eddy-turnover times. Depending on ini- 
tial conditions the energy spectrum at 1024 -resolution 
is still transient at that time. 

In case II, pictured in Fig. [21 strong anisotropy is 
generated due to turbulence depletion along the mean 
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FIG. 2: Field-perpendicular total (solid), kinetic (dashed), 
and magnetic (dotted) energy spectra (normalized, time- 
averaged, and compensated) in 1024 2 x 256 case II simula- 
tion with bo = 5. Dash-dotted curve: high-k part of field- 
parallel total energy spectrum. Inset: perpendicular total en- 
ergy spectrum for resolutions of 512 2 (dash-dotted) to 1024 2 
(solid). 

magnetic field, b , (cf. also [3 E HJ |U ) . This 
is visible when comparing the normalized and time- 
averaged field-perpendicular one-dimensional spectrum, 
Ek ± = J J dkidk2E(k±, hi, fe) (solid line) with the field- 
parallel spectrum, defined correspondingly and adum- 
brated by the dash-dotted line in Fig. |3 The fixed k±- 
axis is chosen arbitrarily in the fci-fc2-plane perpendicu- 
lar to bo where fluctuations are nearly isotropic. For the 
strong b chosen here, field-parallel and -perpendicular 
energy spectra do not differ notably from the ones found 
by considering the direction of the local magnetic field as 
done e.g. in P^El . The field-parallel dissipation length 
is larger than in field-perpendicular directions because of 
the stiffness of magnetic field lines. The numerical reso- 
lution in the parallel direction can, therefore, be reduced. 

While there is no discernible incrtial range in the par- 
allel spectrum, its perpendicular counterpart exhibits an 
interval with Iroshnikov-Kraichnan scaling, E^ x ~ k , 
(Note that due to identical energy scales in Figs. 
and the absolute difference between Kolmogorov and 
Iroshnikov-Kraichnan scaling is the same as in Fig^l. 
This is in contradiction to the anisotropic cascade phe- 
nomenology of Goldreich and Sridhar for strong turbu- 
lence predicting E^ ± ~ k ± °^ 3 7] and with numerical 
studies claiming to support the GS picture How- 
ever, the strength of bo in these simulations is of the or- 
der of the turbulent fluctuations and consequently much 
weaker than for the anisotropic system considered here. 
We note that indication for field-perpendicular IK scaling 
has been obtained in earlier simulations at lower resolu- 
tion using a high-order hyperviscosity and with a stronger 
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mean component, &o/i~3x 10 2 |25| . The authors of the 
aforementioned paper, however, are unsure whether they 
observe a numerical artefact or physical behavior. 

The strongly disparate spectral extent of field-parallel 
and -perpendicular fluctuations suggests that Alfven 
waves propagating along the mean field do not have a 
significant influence on the perpendicular energy spec- 
trum (in the sense of Goldreich-Sridhar, cf. also |2l|). 
Instead, the strong bo constrains turbulence to quasi-two- 
dimensional field-perpendicular planes as is well known 
and has been shown for this particular system |isj . 

Another intriguing feature of system II is that Ef ~ 
E^f with only slight dominance of E M (cf. Fig. EJ in 
contrast to the growing excess of spectral magnetic en- 
ergy with increasing spatial scale for case I. Since both 
states are dynamically stable against externally imposed 
perturbations (as has been verified numerically), they 
presumably represent equilibria between two compet- 
ing nonlinear processes: field-line deformation by tur- 
bulent motions on the spectrally local time scale tnl ~ 
£/vi ~ (k 3 Ef) leading to magnetic field amplifica- 
tion (turbulent small-scale dynamo) and energy equipar- 
tition by shear Alfven waves with the characteristic time 
t~A ~ £/bo ~ (kbo) _1 (Alfven effect). The conjecture 
can be verified via the EDQNM closure approximation 
|2fij which yields evolution equations for kinetic and mag- 
netic energy spectra |27| by including a phenomenological 
eddy-damping term for third-order moments. The spec- 
tral evolution equation for the signed 32] residual energy, 
E K = E M — E K , in the case of negligible cross helicity 
reads [28j |: 

(d t + (ji + ti) k 2 ) Ef = f A d P d q e kpq (t* + t* + t*$) 

with the spectral energy flux contributions 

r£ = m kpq ^EfEf + r kpq ^EfEf , 
T« s = -m kpqP E q Ef ~ t kpqP EfE k , 
Tgyn = ^(k 2 E p E q ~p 2 E q E k ). 

The geometric coefficients m kpq , r kpq , s kpq , t kpq , a con- 
sequence of the solenoidality constraints , are given in 
|28j . The 'A' restricts integration to wave vectors k, p, 
q which form a triangle, i.e. to a domain in the p-q plane 
which is defined by q = \p + k\. The time Q kpq is char- 
acteristic of the eddy damping of the nonlinear energy 
flux involving wave numbers k, p, and q. It is defined 
phenomenologically but its particular form does not play 
a role in the following arguments. 

Local triad interactions with k ~ p ~ q are dominating 
the hydrodynamic turbulent energy cascade and lead to 
Kolmogorov scaling of the associated spectrum (cf., for 
example, |22|). In contrast, the nonlinear interaction of 
Alfven waves includes non-local triads with, e.g., k <C 
p ~ q. In this case a simplified version of equation 



can be derived: 

d t Ef = -T k kEf ee T« lf , (5) 

with r fe = ffc AqQ kpq Ef H3 ~ kE M Q. 

It is now assumed that the right hand side of iQJ) can 
be written as Tf u + T* yn ftjfr. This states that the resid- 
ual energy is a result of a dynamic equilibrium between 
turbulent dynamo and Alfven effect. For stationary con- 
ditions and in the inertial range, dimensional analysis of 
(@J and © yields k z E\ ~ k 2 E M Ef which can be re- 
written as 

Ef ~ kE 2 k . (6) 

The relaxation time, 0, appears as a factor on both sides 
of the relation and, consequently, drops out. We note 
that with ta ~ (kbo) -1 , where bo is the mean magnetic 
field carried by the largest eddies, bo ~ (E^) 1 / 2 , and by 
re-defining r N L ~ (/(v^+bj) 1 ^ 2 ~ (fc 3 ^) -1 / 2 (for system 
II all involved quantities are based on field-perpendicular 
fluctuations) relation JHl can be obtained in the physi- 
cally more instructive form 

^ R -f— YE,. (7) 

The modification of tnl is motivated by considering that 
gradients of the Alfven speed contribute to nonlinear 
transfer as much as velocity shear (see, e.g., [30jV 

For the examined setups relation @ is consistent with 
the underlying physical idea of dynamical equilibrium be- 
tween Alfven and dynamo effect. At small scales with 
k 3> ko (for system II: ko ~ kf), Alfvenic interaction 
always dominates the energy exchange since ta <C tnl 
(e.g. at k = 0.3/p 1 for system I: ta — 5x 10~ 2 , t N l — 0.2, 
for system II: ta — 1 X 10~ 2 ,tnl — 0.1) which re- 
sults in approximate spectral equipartiton of kinetic and 
magnetic energy. At larger spatial scales the Alfven ef- 
fect becomes less efficient in balancing the transforma- 
tion of kinetic to magnetic energy by the small-scale dy- 
namo with ta — tnl (e.g. at k = O.Ol/^for system 
I: ta ~ 0.9, t nl — 0.8, at k = 3 x lO -3 /^ 1 for system 
II: ta — 1.2, tnl — 0.9) allowing larger deviations from 
equipartition. 

An interesting consequence of ijBJ is that the differ- 
ence between possible spectral scaling exponents, which 
is typically small and hard to measure reliably, is en- 
larged by a factor of two in Ef. Even with the limited 
Reynolds numbers in today's simulations such a magni- 
fied difference is clearly observable (e.g. dash-dotted lines 
in Figs. Q]and|2J). For system I with Kolmogorov scaling, 
E k ~ fc~ 5 / 3 (Fig. HJ, relation © predicts Ef ~ k~ 7 / 3 
in agreement with the simulation (Fig. EJ) . In the case of 

3/2 

Iroshnikov-Kraichnan behavior, E k± ~ kj_ as realized 
in system II (Fig. |2J), Ef ~ kj 2 is obtained. This result 
is confirmed by the residual energy spectrum shown in 



4 



1.000 




FIG. 3: Compensated, space-angle-integrated residual energy 
spectrum, E^, for same system as in Fig. Q Dash-dotted 
line: fc _2 -scaling. 
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FIG. 4: Compensated field-perpendicular residual energy 
spectrum for the same system as in Fig. [5] 

Fig. 0] (cf. also [U for two-dimensional MHD simula- 
tions and for spectral model calculations). 

In summary, based on the structure of the EDQNM 
closure equations for incompressible MHD a model of 
the nonlinear spectral interplay between kinetic and mag- 
netic energy is formulated. Throughout the inertial range 
a quasi-equilibrium of turbulent small-scale dynamo and 
Alfvcn effect leads to the relation, Ef ~ kE\, linking to- 
tal and residual energy spectra, in particular ~ fc -7 / 3 
for E k ~ k- r °/ 3 and Ef ~ fc~ 2 for E k ~ fc~ 3 / 2 . Both 
predictions are confirmed by high-resolution direct nu- 
merical simulations, limiting the possible validity of the 
Goldreich-Sridhar phenomenology to MHD turbulence 
with moderate mean magnetic fields. 

The authors would like to thank Jacques Leorat and 
Dieter Biskamp for helpful discussions. WCM acknowl- 



edges financial support by the CNRS and CIAS, Paris 
Observatory. 



[1] S. Ortolani and D. D. Schnack, Magnetohydrodynamics of 
Plasma Relaxation (World Scientific, Singapore, 1993). 

[2] Y. B. Zeldovich, A. A. Ruzmaikin, and D. D. Sokoloff, 
Magnetic Fields In Astrophysics (Gordon and Breach 
Science Publishers, New York, 1983). 

[3] D. Biskamp, Magnetohydrodynamic Turbulence (Cam- 
bridge University Press, Cambridge, 2003). 

[4] A. N. Kolmogorov, Proceedings of the Royal Society A 
434, 9 (1991), [Dokl. Akad. Nauk SSSR, 30(4), 1941]. 

[5] P. S. Iroshnikov, Soviet Astronomy 7, 566 (1964), [As- 
tern. Zh., 40:742, 1963]. 

[6] R. H. Kraichnan, Physics of Fluids 8, 1385 (1965). 

[7] P. Goldreich and S. Sridhar, Astrophysical Journal 485, 
680 (1997). 

[8] S. Sridhar and P. Goldreich, Astrophysical Journal 432, 
612 (1994). 

[9] W.-C. Miiller and D. Biskamp, in Turbulence and Mag- 
netic Fields in Astrophysics, edited by E. Falgarone and 
T. Passot (Springer Berlin, 2002), vol. 614 of Lecture 
Notes in Physics, pp. 3-27. 

[10] R. Grappin, A. Pouquet, and J. Leorat, Astronomy and 
Astrophysics 126, 51 (1983). 

[11] A. Vincent and M. Meneguzzi, Journal of Fluid Mechan- 
ics 225, 1 (1991). 

[12] D. Biskamp and W.-C. Miiller, Physical Review Letters 

83, 2195 (1999). 

[13] W.-C. Miiller and D. Biskamp, Physical Review Letters 

84, 475 (2000). 

[14] N. E. L. Haugen, A. Brandenburg, and W. Dobler, Phys- 
ical Review E 70, 016308 (2004). 

[15] R. J. Leamon, C. W. Smith, N. F. Ness, W. H. 
Matthaeus, and H. K. Wong, Journal of Geophysical Re- 
search 103, 4775 (1998). 

[16] M. L. Goldstein and D. A. Roberts, Physics of Plasmas 
6, 4154 (1999). 

[17] Y. Kaneda, T. Ishihara, M. Yokokawa, K. Itakura, and 

A. Uno, Physics of Fluids 15, L21 (2003). 
[18] W.-C. Miiller, D. Biskamp, and R. Grappin, Physical 

Review E 67, 066302 (2003). 
[19] R. Grappin, Physics of Fluids 29, 2433 (1986). 
[20] J. V. Shebalin, W. H. Matthaeus, and D. Montgomery, 

Journal of Plasma Physics 29, 525 (1983). 
[21] R. M. Kinney and J. C. McWilliams, Physical Review E 

57, 7111 (1998). 
[22] S. Oughton, W. H. Matthaeus, and S. Ghosh, Physics of 

Plasmas 5, 4235 (1998). 
[23] J. Cho, A. Lazarian, and E. T. Vishniac, Astrophysical 

Journal 564, 291 (2002). 
[24] J. Cho and E. T. Vishniac, Astrophysical Journal 539, 

273 (2000). 

[25] J. Maron and P. Goldreich, Astrophysical Journal 554, 
1175 (2001). 

[26] S. A. Orszag, Journal of Fluid Mechanics 41, 363 (1970). 

[27] A. Pouquet, U. Frisch, and J. Leorat, Journal of Fluid 
Mechanics 77, 321 (1976). 

[28] R. Grappin, U. Frisch, J. Leorat und A. Pouquet, As- 
tronomy and Astrophysics 105, 6 (1982). 



■5 



[29] M. Lesieur, Turbulence in Fluids (Kluwer Academic Pub- 
lishers, Dordrecht, 1997). 

[30] J. Heyvaerts and E. R. Priest, Astronomy & Astrophysics 
117, 220 (1983). 

[31] D. Biskamp, Chaos, Solitons & Fractals 5, 1779 (1995). 

[32] The other definition of involving the modulus op- 



erator avoids case differentiations since the applied di- 
mensional analysis is unable to predict the sign of E R . 
However, the physical picture underlying Eqs. (6) and (7) 
implies E^f > E^ as it expresses an equilibrium between 
magnetic energy amplification and equipartition of E K 



